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Microscopic flows are almost universally linear, laminar and stationary because Reynolds number, 
Re, is usually very small. That impedes mixing in micro-fluidic devices, which sometimes limits their 
performance. Here we show that truly chaotic flow can be generated in a smooth micro-channel of a 
uniform width at arbitrarily low Re, if a small amount of flexible polymers is added to the working 
liquid. The chaotic flow regime is characterized by randomly fluctuating three-dimensional velocity 
field and significant growth of the fiow resistance. Although the size of the polymer molecules 
extended in the flow may become comparable with the micro-channel width, the flow behavior is 
fully compatible with that in a table-top channel in the regime of elastic turbulence. The chaotic 
flow leads to quite efficient mixing, which is almost diffusion independent. For macromolecules, 
mixing time in this microscopic flow can be three to four orders of magnitude shorter than due to 
molecular diffusion. 
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PACS numbers: 47.27.-1,47.50. +d,83.50.-v 
I. INTRODUCTION. 

Flows of liquids in microscopic channels have been at- 
tracting increasing interest recently due to fast develop- 
ment of microfluidics and soft lithography [1,2]. The mi- 
crofluidic devices, which are becoming increasingly ad- 
vanced and reliable, allow dramatic reduction of amounts 
of reagents required for fine chemistry and biochemistry 
[3], well controlled manipulation and sophisticated ex- 
periments on individual cells [4-6] and macromolecules 
[7]. The microscopic flows are almost universally lami- 
nar, with linear dependence of the flow rate on the driv- 
ing force. They also usually remain steady as long as 
the driving force does not change. All of that has to 
do with low to moderate values of the Reynolds number. 
Re = Vdp/rj, which is a general measure of non- linear in- 
ertial effects in the flow and of likelihood to find it being 
chaotic or turbulent [8]. Here V is the flow velocity, d is 
the diameter of the channel, and p and 77 are the density 
and the viscosity of the fluid, respectively. When d is 
reduced, the flow velocity needed to reach a given high 
Re, which is required to generate chaotic or turbulent 
flow, increases as d~^ . The driving pressure per unit 
length scales like AP/AL ^ ijV/d'^ giving AP ^ AL/d^ 
at a given Re in the channel. Therefore, if the chan- 
nel proportions are preserved, AP grows quadratically 
with d^^ and when the channels are only a few tens of 
microns wide, achieving sufficiently high Re requires im- 
practically high driving pressures. 

The laminar character of microscopic flows has multi- 
ple practical advantages including possibility of precise 
control of flow velocity, chemical concentration proflles 
[4-6] , and targeted delivery of chemicals and particles [5] . 



On the other hand, the laminar flows have an inherent 
problem of inefflcient mass transfer in directions perpen- 
dicular to the main flow, which occurs due to molecu- 
lar diffusion only. Diffusion time, d^ /D, across a typical 
micro-channel with a width of 100 fim is on the order 
of 100 s even for moderate size proteins, such as bovine 
serum albumin with a diffusion coefficient of D w 3 • 10~^ 
cm^/s in water [9]. 

A few techniques have been proposed to generate stir- 
ring by a three-dimensional flow in order to increase the 
rate of mixing in the micro-channels. They include ap- 
plication of time dependent external forces [10,11] and 
raising Re to moderately high values in curvilinear three- 
dimensional channels [12,13]. An ingenious method of 
mixing has been suggested recently, which involves spe- 
cial "herring-bone" patterning of a micro-channel wall 
to generate fluid motion perpendicular to the main flow 
direction in the linear, low Re regime [14]. The fluid 
elements are continuously stretched and folded in the 
flow as they advance along the channel. That separates 
closely spaced fluid particles and brings distant particles 
together, dramatically reducing the characteristic length 
scales and diffusion time and increasing the rate of ho- 
mogenization of the mixture. The flow was stationary 
in the laboratory frame, however. Therefore, the con- 
centration profiles of the fluorescent dye used as a tracer 
were uniquely defined by the entrance conditions and the 
channel geometry [14] . They did not change in time, and 
there was always some constant difference in concentra- 
tion between neighboring points. 

The basic condition of linearity in low Re flows can be 
changed by adding flexible high molecular weight poly- 
mers to the working hquid [15]. Solutions of those poly- 
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mers arc known as non-Newtonian visco-clastic fluids 
[15]. Mechanical stress in these fluids depends on the 
flow history with some characteristic relaxation time, A, 
which for dilute solutions is a time of relaxation of in- 
dividual polymer molecules. Another specific property 
of the polymer solutions is the non-linear dependence of 
the polymer contribution to the stress on the rate of de- 
formation in the flow, W [15]. This non-linearity usu- 
ally becomes significant, when the Weissenberg number, 
Wi = XVV, becomes on the order of unity. The non- 
linear growth of the elastic polymer stresses is especially 
striking in extensional flows at Wi > 1/2, where appar- 
ent viscosity of polymer solutions can rise by up to three 
orders of magnitude as the total deformation increases 
[16]. In pure shear flows the major non-linear elastic ef- 
fect shows up in the appearance of negative normal stress 
along the flow direction. It leads to the well known effect 
of rod climbing [15,17] and causes purely elastic instabil- 
ities in curvilinear flows of viscous polymer solutions in 
table-top set-ups, when inertial effects are virtually ab- 
sent [18]. The non-linear polymer stresses in curvilinear 
shear flows can also lead to elastic turbulence, a random 
multi-scale three-dimensional flow, which can develop at 
arbitrarily low Re [19]. Elastic turbulence causes sharp 
growth of the flow resistance [19,20], and it was foimd to 
generate efficient mixing in a table-top curvilinear chan- 
nel [21]. 

Recently it was shown that the non-linear elasticity of 
the polymer solutions can cause a purely elastic transi- 
tion and non-linear growth of flow resistance in a mi- 
croscopic channel with contractions [22]. Those non- 
linear effects, however, were due to regions of fast ex- 
tensional flow near the contractions, where individual 
polymer molecules are expected to partially unravel at 
Wi > 1/2 [23,24]. The non-linearity in resistance was 
the most dramatic feature of the transition and it became 
quite signiflcant before any substantial changes in the 
flow pattern could be seen [22] . Therefore the non- linear 
resistance could be understood as a simple additive effect 
of individual molecules forced through the contractions. 
Although the flow became rather irregular at higher flow 
rates, the fluctuating flow regions were mostly near the 
contractions, and no detailed study of those fluctuations 
was made. 

If the basic linear flow is a pure shear as in Ref. [19-21] , 
a non-linear elastic transition can only occur through 
a major reorganization of the flow structure. In Ref. 
[19-21] the secondary flow generated above the instabil- 
ity threshold had well expressed turbulent features and 
involved irregular fluid motion in a broad range of tem- 
poral and spatial scales. That implies an essentially col- 
lective effect of the polymer molecules on the flow, and 
the polymer solution behaving as a visco-elastic contin- 
uum. If the size of the set-up is reduced to a micro-scale, 
non-homogeneity in the polymer concentration and ex- 
tension of polymer molecules stretched in the flow may 



become comparable with the size of the set-up, and larger 
than size of some of the generated vortices. Therefore, 
whether or not the purely elastic instability and elastic 
turbulence can be reproduced in a microscopic shear flow 
is still an open question. 

Here we show that a fully developed chaotic flow sim- 
ilar to the clastic turbulence can be generated in a flat 
curvilinear microscopic channel with smooth walls and 
uniform width at arbitrarily low Re, if the working liq- 
uid contains a small amount of high molecular weight 
polymers. The flow is characterized by significant non- 
linear growth in resistance, randomly fluctuating velocity 
field and chaotic three-dimensional mixing patterns. It is 
further shown that stirring by the flow results in efficient 
mixing in the micro-channel with characteristic mixing 
length significantly shorter compared with the "herring- 
bone" patterning method reported before [14]. The char- 
acteristic mixing times for solutions of macromolecules 
are reduced by three to four orders of magnitude com- 
pared with molecular diffusion. 

II. MATERIALS AND METHODS. 

A. Device fabrication. 

The micro-channel devices consist of a silicon elas- 
tomer (Sylgard 184 by Dow Corning) chip sealed to a #1 
microscope cover glass. The channel structure of the chip 
was fabricated using the technique of soft lithography 
[1]. First, a negative master mold was fabricated in UV- 
curable epoxy (SOTEC micro-systems SU8-1070) by us- 
ing conventional photolithography. The epoxy was spun 
onto a silicon wafer at 1800 rpm for 60 s to create a 100 
jum layer and patterned by using a high-resolution nega- 
tive transparency mask. Liquid elastomer was poured on 
the mold to a thickness of « 5 mm and cured in an 80C 
oven for 1 hr 30 min. After that the elastomer was peeled 
off the mold, trimmed to its final size and liquid feeding 
ports were punched by using a 20-gauge lucr stub. The 
patterned side of the chip was bonded to the cover glass 
by overnight baking in the 80C oven. 

B. Experimental set-up. 

A snapshot of our first microfiuidic set-up is shown in 
Fig.l. It has a uniform thickness d = lOO^m. Its main 
active element is a curvilinear channel with square cross- 
section. It is a chain of 40 identical segments, which 
are couples of interconnected half-rings with inner and 
outer radii Ri = 100/Lxm and Ro = 200/[xm, respectively. 
Fig. lb. It has the same proportions as the table- top 
channel, which was used in the elasticity induced mix- 
ing experiments reported before [21], but its dimensions 
are reduced by a factor of 30. Because of the periodic 
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structure of the channel, it is convenient to use the num- 
ber of a segment, N, starting from the inlet as a discrete 
linear coordinate along the channel. The auxiliary recti- 
linear channel (b) has width of 90 /im and total length of 
about 72.5 mm. Channel (b) and the comparator region 
(c) serve to make differential in situ measurements of flux 
vs. pressure by the method described in Ref. [22]. 

C. Flow control. 

The flow in the micro-channels was generated and con- 
trolled by pressure differences between the inlets and the 
outlet. Fig. la. The pressures were generated hydrostati- 
cally using long vertical rails with precise rulers and slid- 
ing stages. Working liquids were kept in 30 ml plastic 
syringes, which were held upright, open to the atmo- 
sphere and connected to the two inlets and the outlet by 
plastic tubing with internal diameter of 0.76 mm. The 
pressure drop in the tubing was estimated below 1% of 
the total. The two syringes feeding the inlets were at- 
tached to the sliding stages. The difference in liquid ele- 
vation between these two syringes and the outlet syringe 
was measured and adjusted with a precision of about 
0.1 mm corresponding to 1 Pa in pressure. Dependence 
of the volumetric low rate, Q, in the curvilinear chan- 
nel on the pressure difference between inlet 1 and outlet, 
Fig.l, was determined with a relative precision of about 
0.5% using an in-situ measurement technique described 
elsewhere [22]. A syringe pump (PHD 2000 by Harvard 
Apparatus Inc.) with a 50 ml gastight Hamilton syringe 
was used for an absolute flow rate calibration. 



D. Polymer solutions. 

The polymer used was polyacrylamide, PAAm, by 
Polysciences Inc. with high molecular weight Mu, = 
1.8 • 10^, which was the same polymer sample as in Rcf. 
[21]. It was dissolved at identical concentrations of 80 
ppm by weight in two Newtonian solvents with different 
viscosities. The solvent for the low viscosity solution 1 
was a 35% solution of sucrose in water with 1% of NaCl 
added to fix the ionic contents. The Newtonian viscosity 
of the solvent, ry^, was 4.2 mPas at the room temperature 
of 22C. The viscosity of solution 1, 77, was 5.6 mPas at 
a shear rate of 50 s~^, suggesting a dilute polymer solu- 
tion. The Newtonian solvent for the high viscosity solu- 
tion 2 was sugar syrup containing 64.4% sucrose and 1% 
of NaCl in water. The viscosity of solvent 2 was 114 mPas 
at 22c, and the viscosity of solution 2 was 138 mPas at 
a shear rate of 2 s~^. The polymer relaxation time, A, 
of solution 2 measured by phase shift between stress and 
shear rate in an oscillatory flow regime was 1.1 s. The 
measurements of viscosity and relaxation time were made 



using a high precision rhcometer (ARIOOO by TA Instru- 
ments). Relaxation time of solution 1 was estimated as 
0.04 s with the assumption that A scales linearly with r]s 
[15]. 

The overlap concentration c*, taken as a concentration 
at which the viscosity ratio reached r]/r]s = 2, was 200 
ppm by weight for solvent 2, corresponding to molec- 
ular concentration of n = 8.76 • 10^^ cm~"^. Charac- 
teristic size of the polymer coils at rest can be esti- 
mated from this as n"^/"^ w 0.5/im, and characteristic 
distance between them at 80 ppm by weight can be esti- 
mated as 0.7 /xm. These estimates are well supported by 
the data on PAAm taken from the literature [25]. One 
can use the Mark-Houwink scaling relation for PAAm, 
[ri] = 6.31 X 10-^M°-^ (in ml/g) [25], where [r]] is the in- 
trinsic viscosity of the polymer defined as [77] — [c^o- 
At molecular weight Myj = 1.8 x 10'' one gets [rj] = 4020 
ml/g. Defining c* = l/[ri] [15] we obtain c* = 250 ppm 
by weight for an aqueous solution in a good agreement 
with the above estimate. Further, we can compare the 
estimated size of the PAAm coils with the data obtained 
from light scattering Rcf [25]. Plugging in M^^ = 1.8x10'' 
into an interpolation relation from Ref. [15], one gets 
Rg w 0.4/im for the radius of gyration of the coils, which 
is rather consistent with the above estimate obtained 
from c*. 

It is worth noting that c* = 200 ppm is very close to 
the value found for A-phage DNA [26] , which has a com- 
parable molecular weight of 3.1 x 10''. Radius of gyration 
of the A-phage DNA was found to be 0.73/xm [27], which 
is quite in line with the above estimate for the PAAm coil 
size. The full contour length of a PAAm molecule having 
= 1.8 X 10'' and thus consisting of 2.5 x 10^ monomers 
(having molecular weight 71.08 g/mol) can be estimated 
as about 50/im, if a monomer length of 0.2 nm is as- 
sumed. It is significantly larger than the contour length 
of the A-phage DNA (which is equal to about 16/im) and 
twice smaller than the micro-channel diameter. 

The experiments on mixing were carried out with poly- 
mer solution 2, and then fluorescent dyes with low dif- 
fusivities were added to the solution and used as trac- 
ers. Those were a few different samples of fluorescein- 
conjugated Dextran, FITCD, by Sigma with average 
molecular weights, M, varying from 10 kDa to 2 MDa. 
In spite of the relatively high molecular weight of FITCD 
it did not have any measurable influence on the solu- 
tion rheology due to high rigidity of the polysaccharide 
molecules. Diffusion coefflcients of the FITCD samples 
in water were estimated using the data in Ref. [28] , giving 
values from 9.1-10"'' to 7.4-10"^ cm /s that corresponded 
to a broad range of biological macromolecules. The dif- 
fusion coefficients in solvent 2 were estimated with as- 
sumption that D ^ I/t/s, resulting in Di = 6.6- 10^^ and 
D2 = 5.4 • 10-1° cmVs for 10 kDa and 2 MDa FITCD, 
respectively. 
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E. Measuring flow velocity. 

Measurements of the flow velocity in the micro-channel 
were carried out using custom developed microscopic par- 
ticle image velocimctry, micro-PIV. The polymer solution 
was seeded with 0.2 fim yellow-green fluorescent beads 
(Polysciences) , and epi- fluorescent imaging of the flow 
in the micro-channel, Fig. 1, was made with an inverted 
microscope (Olympus IMT2) and narrow band excitation 
and emission filters in the dichroic filter cube. The objec- 
tive was a LWD 20 x, N.A. = 0.40 , and the images were 
projected onto a CCD array with 640 x 480 pixel resolu- 
tion (PixelFly camera by PCO, Germany) and digitized 
to 12 bits. The snapshots of the flow were taken with even 
time intervals of 40 ms, and digitally post-processed. Im- 
ages of out-of-focus particles were disregarded. Velocity 
field was found by cross-correlating positions of the parti- 
cles in two consecutive snapshots, and the particle veloc- 
ity vectors were neighbor-validated. (The calculated ve- 
locity field corresponded to the time interval between the 
two snapshots.) The collected time series represented ve- 
locity values measured at equal distances from intercon- 
nections of two half-rings a.t N = 35, and averaged over 
a 20 X 20/Ltm square region at the middle of the channel 
and over 4 /xm across the channel mid-plane. 

F. Measuring tracer concentration profiles. 

Concentrations of the fluorescent dyes, which were 
used as passive tracers in the experiments on mixing in 
the channel, were measured using a commercial confo- 
cal microscope (Fluoview FV500 by Olympus). It was 
equipped with a 40x N.A.= 0.85 infinity corrected ob- 
jective and a 12 bit photomultiplicr. The scanning was 
done at a rate of 56 lines per second and 512 pixels per 
line corresponding to a step of 0.18 /im per pixel . 

III. RESULTS. 

We measured volumetric flux rate, Q, of solution 1 
through the curvilinear channel, Fig.l, in a broad range 
of applied pressures and calculated the resistance factor, 
Z — P/Q, where AP is pressure drop per segment. The 
resistance factor is a constant proportional to viscosity 
for Newtonian fluids in linear, low Re regime, and it can 
be used as a measure of turbulent flow resistance in large 
channels at high Re. Fig. 2 shows dependence of Z on 
Q for solution 1, after Z is divided by a resistance factor, 
Zo, found for a Newtonian liquid with the same viscosity, 
rj. The ratio Z/Zq is constant equal to unity in the linear 
regime at low Q. At Q of about 8.5 nl/s, however, a non- 
linear transition occurs; Z/Zq starts to grow and reaches 
a factor of about 2.8 at high Q. The Reynolds number 



for the channel flow can be dcflncd as Re = Qp/{rjd) . 
It was 0.017 at the transition point and 0.14 at highest 
Q, that we tried, so that the inertial effects were always 
negligible. The Weissenbcrg number can be dcflncd as 
Wi = 4:XQ/(P, and its estimated value at the transition 
point is 1.4, which is comparable with Wi found at purely 
elastic transitions in macroscopic set-ups [18,19,21,29]. 

In order to get detailed information about structure of 
the flow above the non-linear transition we used solution 
2 with high viscosity and large A. The polymer relax- 
ation time deflnes both characteristic time of changes in 
the flow and the inverse of flow velocity in the elastic 
non-linear regime [15,29]. Therefore, the non- linear flow 
in solution 2 was expected to be much slower, and mea- 
surements of its characteristics were expected to be more 
feasible with the standard video microscopy techniques. 
Using the micro-PIV we measured flow velocity in the 
middle of the curvilinear channel. Dependence of RMS 
of fluctuations of the longitudinal component of the flow 
velocity, V{"^'', (which is the velocity component along 
the main flow) on AP is shown in Fig. 3. The fluctuations 
are virtually absent in the linear regime at low pressure. 
At AP of about 50 Pa, however, Vi starts to fluctuate, 
and V{^^ begins to grow quickly and non-linearly. It 
can be learned from the inset in Fig. 3 that the longitu- 
dinal component of the flow velocity, Vi , grows linearly 
at low AP, but its growth slows down at the same crit- 
ical AP of about 50 Pa, which is another evidence for 
a non-linear elastic flow transition taking place. At the 
transition point, the average longitudinal flow velocity, 
estimated from the micro-PIV measurements in different 
points across the channel, was Vi w 80/im/s. That gives 
estimates of Q « 0.8 nl/s, Pe « 8 • 10"^ and Wi w 3.5 
for the elastic non-linear transition. This value of the 
Weissenberg number is very close to Wic ~ 3.2 found for 
the transition to chaotic flow in the table-top curvilinear 
channel in Ref. [21] 

A typical time series of Vi above the transition is shown 
in Fig. 4a. The velocity is strongly fluctuating and its 
time dependence has well expressed chaotic appearance. 
The chaotic character of the velocity fluctuations is con- 
flrmed by analysis of its time correlations. The velocity 
autocorrelation function shown in Fig. 4b does not have 
distinct peaks and decays uniformly. 

For the experiments on mixing the design of the micro- 
channel was slightly modifled to enable side by side in- 
jection of two streams of solution 2, one with and one 
without FITCD, to the channel inlet. Fig. 5a. Apart from 
Co = 280 ppm by weight of FITCD added to one of them 
the polymer solutions were identical and were injected 
at equal flow rates by careful adjustment of the driving 
pressures, Fig. 5a. The set-up was flrst tested with the 
plain solvents without PAAm added. The flow appeared 
laminar at all AP that we applied, and the interface be- 
tween the streams with and without FITCD remained 
smooth and sharp along the whole channel with only mi- 
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nor smearing by diffusion, Fig. 5b. 

The situation was similar with the polymer solutions 
in the linear regime at low AP. However, when the driv- 
ing pressure was raised above the non-linear transition 
threshold, fluctuating flow velocity produced significant 
stirring and complex and chaotically changing tracer con- 
centration profiles, Fig. 5c, d. We studied mixing in 
the channel in detail at AP = 134 Pa, corresponding 
to a flow rate about twice above the non-linear transi- 
tion threshold (cf. Rcf. [21]) and Vi of about 173 /im/s. 
Variation of the tracer concentration profiles with time 
at different distances from the inlet is illustrated by the 
space-time plots in Fig. 6a, b. One can observe that the 
tracer concentration appears to fluctuate quite randomly 
without any apparent scale in time or space. Next, one 
can see in Fig. 6a, taken at iV = 12, that the left side of 
the channel, where the tracer was initially injected, looks 
much brighter and has much higher average concentra- 
tion of the tracer. Although also noticeable in Fig. 6b 
taken further downstream, at iV = 18, this feature is 
clearly weaker there. Thus, stirring by the fluctuating 
velocity field seems to create a more symmetric distri- 
bution of the tracer between the two sides of the chan- 
nel. In order to validate this observation, we measured 
time averages of the tracer concentration, c, at different 
positions across the channel and at different N, Fig. 7a. 
One can see that the cross-channel distribution of c/cq 
close to the inlet, at TV = 7, is strongly influenced by 
the asymmetric conditions at the channel entrance. As 
one can learn from the curve at N = 11, however, the 
imprint of the initial conditions is clearly fading as the 
liquid advances downstream and being stirred. Further 
downstream, at 41, asymmetry in the tracer dis- 

tribution introduced by the initial conditions disappears 
completely. Fading of the initial condition influence with 
time and restoration of symmetry in flow in statistical 
sense are both distinct features of chaotic and turbulent 
flows. Therefore, the curves in Fig. 7a provide further 
evidence for truly chaotic nature of the flow in the micro- 
channel. 

A natural parameter characterizing inliomogeneity of 
the mixture is standard deviation of the instantaneous 
local tracer concentration, c, from its overall average 
value (c) = cq. It is convenient to divide it by (c) and 
to introduce a dimensionless standard deviation Cstd = 
\J {{c — (c))^) / (c) . At the channel entrance Cstd is equal 
to unity, and it becomes zero, when the liquid is perfectly 
mixed and homogeneous. Dependence of Cstd on N for 
two tracers with diffusion coefficients of Di = 6.6 • 10~^ 
and D2 = 5.4 • 10~^° cm2/s is shown in Fig. 7b in semi- 
logarithmic scale. One can see that the mixture becomes 
increasingly homogeneous as the liquid advances down- 
stream, and the parameter Cstd decays exponentially with 
A'' for the both tracers. An exponential decay was also 
found in the case of elastic turbulence in a macro-channel 
[21], and it agrees very well with theoretical predictions 



for the so-called Batchelor regime [30 33] of mixing. The 
latter corresponds to a fiow, which is chaotic in time but 
essentially "smooth" in space, in the sense that small ed- 
dies are rare, and the main contribution to mixing comes 
from the largest eddies having the size of the whole sys- 
tem. 

The rates of the exponential decay for the two depen- 
dencies shown in Fig. 7b were 0.217 and 0.137, corre- 
sponding to AN of 4.61 and 7.30, and lengths Li = 4.34 
mm and L2 = 6.88 mm, respectively. The latter can be 
considered as characteristic mixing distances along the 
channel for the two tracers. Characteristic mixing times, 
which can be estimated as tmix = L/V, are then found 
to be tmix,i = 25 s and tmix,2 = 40 s, respectively. In 
the absence of active stirring, mixing would only occur 
through molecular diffusion, with characteristic diffusion 
times across the channel given by tdiff = (P/D resulting 
in tdiff = 1.5-10* s and tdiff, 2 = 1.9-10^ s, respectively. 
Therefore, the stirring produced by the chaotic flow in the 
channel reduces the mixing times for the FITCD macro- 
molecules by three to four orders of magnitude. 

IV. DISCUSSION. 

We studied flow of two dilute polymer solutions in 
a microscopic curvilinear channel. We observed a non- 
linear transition to occur in the flow at Re k, 1.7 • 10^^, 
Wi « 1.4, and Pe « 8 • 10-^ Wi « 3.5, for solution 1 
and solution 2, respectively. The very small and different 
values of Re, and comparable, order of unity values of Wi 
at the transition threshold for the two solutions indicate 
that the transition is of purely elastic nature and can oc- 
cur at arbitrarily low Re. A possible explanation of the 
difference in estimated critical Wi for the two polymer 
solutions is that the actual polymer relaxation time of 
solution 1 may be significantly higher than the estimate 
based on the assumption that A ~ r/g. 
A major global feature of the flow above the non-linear 
elastic transition is fast growth of the resistance. Fig. 2. It 
increases by up to a factor of 2.8 above the resistance for 
a Newtonian fluid with the same shear flow viscosity, 77, 
at the same Q. Since Re is very low, the whole growth of 
the flow resistance should be due to increase in the elastic 
polymer stresses [20] . In a pure shear flow with 7 = 240 
s"'^, corresponding to Q = 60 nl/s (and Wi w 10), viscos- 
ity ratio for polymer solution 1 was rj/rjs ~ 1.22. Hence, 
the average increase in the polymer shear stresses due to 
the secondary flow at Q = 60 nl/s (Fig. 2) can be esti- 
mated as a factor of about 11. Suggesting linear defor- 
mation of the flexible polymer molecules in a shear flow 
at low to moderate Wi [34] , we can estimate extension of 
the molecules at Wi = 10 as 10 • 0.5/im = 5^m. (Flexible 
polymer molecules were found to extend linearly with the 
shear rate, until the extension reached 15-20% of their full 
length [34].) Additional extension of the polymers due to 
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the irregular secondary flow, as suggested by the stress 
growth by the factor of 11, can be estimated as a factor 
of \/TT (suggesting that the secondary flow results in an 
isotropic polymer unravelling and the stress grows as a 
square of the polymer extension [15]). That brings char- 
acteristic size of the extended polymer molecules to 20/um 
range, less than an order of magnitude smaller than the 
diameter of the channel, d = lOO/Ltm. 

The growth of the elastic stresses due to the fluctuat- 
ing secondary flow can only occur as a result of signifi- 
cant reorganization of the flow structure and spontaneous 
generation of regions with strong extensional flow, as in 
the case of the elastic turbulence [19,20]. (In shear flows 
polymer contribution to resistance increases slower than 
linearly with the shear rate, which is called shear thin- 
ning.) This suggestion is quite corroborated by the direct 
flow velocity measurements in solution 2 above the non- 
linear transition. Velocity is found to be strongly fluc- 
tuating. Figs. 3, 4, with RMS of the fluctuations reach- 
ing as much as 10% of the mean longitudinal velocity 
at the center of the channel. Fig. 3, just as in the table- 
top channel in Ref. [21]. The velocity appears to vary 
randomly in time. Fig. 4a, and its auto-correlation func- 
tion decays quite quickly and does not have any distinct 
peaks. Fig. 4b. These all are clear indications of chaotic 
nature of the flow in the channel in the non-linear regime 
at high Wi. Another evidence of chaos comes from the 
experiments on mixing. Figs. 6, 7. Asymmetry in distri- 
bution of c/cq across the channel, which is imposed by 
the conditions at the channel entrance, decays with the 
distance from the entrance, Fig7a. Further, Cstd expo- 
nentially decays with N, as it is supposed to be in the 
chaotic Batchelor regime of mixing [21,32,33]. 

Experimental results presented in Figs. 2-4 and Figs. 6, 7 
and discussed above are fully consistent with a suggestion 
that the regime of elastic turbulence [19] is being real- 
ized in the micro-channel. In fact, the results in Figs. 2-4 
and Figs. 6, 7 agree rather well with the mcasurcincnts in 
macroscopic systems reported before [19,20]. We find 
it rather remarkable that although the size of c^xtendcd 
molecules may become comparable with the channel di- 
ameter, it does not seem to cause any significant new ef- 
fects in the polymer solution dynamics. An essential fea- 
ture of turbulent flows is a broad range of spatial scales, 
at which fluid motion is excited. Unfortunately, limited 
resolution of the micro-PIV technique did not allow us to 
explore properties of the flow down to sufficiently small 
spatial scales and to obtain conclusive data about its spa- 
tial structure. Therefore we can only refer to the flow in 
the micro-channel as being chaotic. 

In the experiments with passive tracers, we found that 
stirring by the chaotic flow results in efficient mixing in 
the channel. Figs. 5-7. So, mixing times for the macro- 
molecules, which we used as tracers, were reduced by 
three to four orders of magnitude compared with pas- 
sive molecular diffusion. Fig. 7b. It is worth noting that 



the characteristic distance of AA^ = 7.3 required for 
mixing of low diffusivity FITCD is comparable with 
A A'' Ri 15 found in the chaotic flow in the table-top chan- 
nel Ref. [21]. We further notice that whereas the diffusion 
time scaled as 1/D, the time of mixing in the chaotic 
flow depended on D very weakly. It increased by only 
60% between the lower and upper curve in Fig. 7b (from 
tmix,i = 25 s and tmix,2 = 40 s), whereas D dropped by 
a factor of 12 (from Di = 6.6 • 10"^ to £>2 = 5.4 • 10"^° 
cm2/s). The weak dependence of tmix on D is quite con- 
sistent with theoretical predictions for chaotic flows and 
the Batchelor regime [30-33]. It implies that the elastic- 
ity induced chaotic stirring in the micro-channel can be 
efficiently used for mixing of liquids with additives of any 
low diffusivity, such as large molecules of DNA, viruses, 
particles, and possibly living cells. 

The relatively long times of mixing in the chaotic fiow 
of solution 2 are due to high viscosity of the solvent, large 
A and low flow velocity at the elastic instability threshold. 
In a much more practically relevant case of water based 
solutions with the viscosity on the order of 0.001 Pas, A 
is expected to be about hundred times lower and the flow 
velocities in the chaotic regime shoiild be about hundred 
times higher [15,29]. Suggesting qualitatively similar flow 
conditions [29] , we expect characteristic mixing times for 
the water based solutions to be on 100 ms scale. Since dif- 
fusion coefficients for macromolecules are proportional to 
l/r]s [15], the ratio between t^iff and tmix for the chaotic 
flow in the water based solutions should be in the same 
range of 10^ to 10^. 

The characteristic mixing length along the channel for 
water based solutions should also remain in the same 
range of about 4-7 mm. It is about 2-3 times shorter 
than characteristic lengths in the "herring-bone" pat- 
terned channel of Ref. [14] at comparable values of Peclet 
number, Pe = Vd/D. Although the method of mixing 
by elasticity induced chaotic flow necessarily requires 
addition of polymers to the working liquid, it may be 
quite practical for many biochemical assays, taking into 
account the very low concentration of the high molecular 
weight polymers used. It does not rely on any special 
patterning of the micro-channels and should be readily 
compatible with rounded channel profiles used for inte- 
grated micro- valves and peristaltic pumps [35], allowing 
efficient mixing in closed loop microscopic flows [36]. 
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FIG. 1. (A) Photograph of the micro-fluidic device. The 
micro-channels were filled with ink for better contrast. (B) 
Magnified image of a section of the functional curvilinear 
channel. 

FIG. 2. Dependence of normalized resistance, Z/Zo, in flow 
of solution 1 through the curvilinear channel on the volumet- 
ric flow rate, Q, in semi-logarithmic scale. 

FIG. 3. Dependence of RMS of fluctuations of the longitu- 
dinal component of flow velocity, Vl"^" , on pressure drop per 
segment, AP, in the center of the micro-channel for solution 
2. Inset: Dependence of time average of the longitudinal flow 
velocity, Vi , in the center of the micro-channel on AP. 

FIG. 4. (A) Time series of the longitudinal flow velocity, 
Vi , in the center of the micro-channel at AP = lOOPa. (B) 
Auto-correlation function for Vi based on about 6000 individ- 
ual velocity measurements. 

FIG. 5. (A) Epi-fluorescent micro-photograph of the en- 
trance area of a micro-channel used in experiments on mix- 
ing. Wide triangular region in front of a curvilinear channel 
allows to adjust equal flow rates for polymer solutions with 
(from below) and without FITCD. (B) Confocal photograph 
of flow in the micro-channel without polymers added. Right 
wall of the channel is shown by a dotted line from below. (C) 
Confocal image of mixing in chaotic flow in the micro-channel 
with solution 2 in it. (D) Confocal imago of cross-section of 
the micro-channel with chaotic flow in solution 2. 

FIG. 6. Space-time plots of FITCD distribution across the 

channel taken at (A) N = 12 and (B) A'" = 18. Confocal scan- 
ning was done along the same line across the channel in the 
mid-plane at equal distances from half-ring interconnections, 
with even time intervals of 0.0177 sec. Profiles of FITCD con- 
centration in consecutive moments of time are plotted from 
top to bottom. 
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FIG. 7. (A) Time average of FITCD concentration, c , di- 
vided by Co, as a function of position, x, along a line across 

the channel taken at different distances from the inlet: N = 7, 
N = 11, and N — 41. The lines across the channel were in the 
mid-plane at equal distances from half-ring interconnections, 
just as for the space-time plots in Fig.6. (B) Standard de- 
viation of FITCD concentration from its average value, Cstd, 
as a function of distance, A'^, from the inlet for two kinds of 
FITCD with average molecular weights of 10 kDa (triangles) 
and 2 MDa (squares). 
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